function plotYieldCurveDecom(yieldDecom,setup)

regime_2  = real(setup.econAct<0);

figure('Name','Term Premia','NumberTitle','off');
dates = datenum(setup.timeIndex);
plotMatYear = [3 5 7 10];
for i=1:size(plotMatYear,2)
    for k=1:size(yieldDecom.matSelect,2)
        if yieldDecom.matSelect(1,k) == plotMatYear(1,i)*12
            subplot(2,2,i)    
            hold on
            timeIndex = str2num(datestr(dates,10))+str2num(datestr(dates,5))/12;
            seriePlot = yieldDecom.termPremia(k,:);
            minY      = min(seriePlot);
            maxY      = max(seriePlot);
            X  = [timeIndex(1,1):1/12:timeIndex(end,1),fliplr(timeIndex(1,1):1/12:timeIndex(end,1))];
            % For shading above zero
            Y1 = [maxY*regime_2' 0*regime_2'];
            fill(X,Y1,[0.8 0.8 0.8]);
            % For shading below zero
            Y2 = [minY*regime_2' 0*regime_2'];
            fill(X,Y2,[0.8 0.8 0.8]);
            
            % plotting the series
            plot(timeIndex,seriePlot,'-k')
            
            ylabel('Basis points')
            hold off
            axis tight
            title(['Term premia: ',num2str(plotMatYear(1,i)), '-year horizon']);   
        end
    end
end

figure('Name','Expected policy rate','NumberTitle','off');
dates = datenum(setup.timeIndex);
plotMatYear = [0.5 1 3 5];
for i=1:size(plotMatYear,2)
    for k=1:size(yieldDecom.matSelect,2)
        if yieldDecom.matSelect(1,k) == plotMatYear(1,i)*12
            subplot(2,2,i)    
            hold on
            timeIndex = str2num(datestr(dates,10))+str2num(datestr(dates,5))/12;
            seriePlot = yieldDecom.rExp(k,:)*100;
            minY      = min(seriePlot);
            maxY      = max(seriePlot);
            X         = [timeIndex(1,1):1/12:timeIndex(end,1),fliplr(timeIndex(1,1):1/12:timeIndex(end,1))];
            % For shading above zero
            Y1 = [maxY*regime_2' 0*regime_2'];
            fill(X,Y1,[0.8 0.8 0.8]);
            % For shading below zero
            if minY < 0
                Y2 = [minY*regime_2' 0*regime_2'];
                fill(X,Y2,[0.8 0.8 0.8]);
            end
            
            % plotting the series
            plot(timeIndex,seriePlot,'-k')
            ylabel('Percentage')
            hold off
            axis tight
            title(['Expected policy rate: ',num2str(plotMatYear(1,i)), '-year horizon']);   
        end
    end
end

end

